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Abstract The CP 2 model, with and without a generalized Hopf term, 
is studied using the collective coordinate approximation. In the spirit of 
this approximation, an ansatz is given which in previous numerical studies 
was seen to give a good parameterization of the numerical solution. The 
equations of motion for the collective coordinates are then solved analytically, 
for solitons close together and for solitons far apart. The solutions show how 
the generalized Hopf term changes the scattering angle which in its absence 
is 90°. 
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1 Introduction 



ex-models in low dimensions have become an increasingly important area of 
research. In two Euclidean dimensions they appear to be the low dimen- 
sional analogues of four-dimensional Yang-Mills theories. They often arise as 
approximate models, in the context of both particle and solid state physics. 
They have recently been used in the construction of models of high T c su- 
perconductivity [HI and of the quantum Hall effect 0. Moreover, they are 
simple examples of harmonic maps studied by differential geometers and, as 
such, are interesting in themselves. In addition, because of the nonlinearity 
of these models, their classical solutions represent structures which resemble 
solitons of (1+1) dimensional models and could become very useful in the 
description of many physical phenomena. 

Recently, interesting scattering processes of these soliton-like objects were 
studied with the help of the collective coordinates approximation 011110 
and with the help of numerical simulations. The numerical simulations were 
performed for the CP 1 model @0, for its modification by the addition of 
potential-like and Skyrme-like terms |§||, and most recently for the CP 2 
model with and without a generalized Hopf term []10|1||11|| . In Ref. 11, also 
the Ginibre et al. existence proof JT2| and the Cauchy-Kowalewskyi theorem 
were used to study CP 2 soliton scattering. Although there is no mathemati- 
cally rigorous justification for the collective coordinate approximation in the 
CP n model yet, Manton's arguments, first put forward for the monopole 
motion |T3L should also apply to the CP 2 soliton motion. In this paper we 



therefore use this approximation to study some aspects of the CP 2 soliton 
scattering. We again find the 90° scattering without, and a deviation from 
the 90° scattering, with the generalized Hopf term. 

The paper is organized as follows. In section 2, we describe the CP 2 
model without and with generalized Hopf term. We also formulate a Cauchy 
problem suitable for the description of solitons which merge at some stage of 
the scattering process. In section 3, the description of the scattering process 
in terms of collective coordinates is given. In section 4, we solve the equations 
for the collective coordinates for solitons close enough together, i.e. for short 
enough times before and after the collision. Solitons far apart are studied in 
section 5. 
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2 The model 



The Lagrangian for the CP 2 model is, 

C = (D"Z a )*(D li Z a ), ^ = 0,1, 2; a = 1,2, 3; (2.1) 

where Z a is a complex function of {t^x 1 ^ 2 ) on M 3 (or a function of the 
variables (t,z = \{x l + ia; 2 ),z*)), and D^Z a = d^Z a - Z* b {d^Z b )Z a . Also, 
the Z a have to satisfy the condition Z* a Z a = 1. a,b,.. are raised and 
lowered with the metric t] = diag(+l, +1, +1). The space-time metric is 
g = dia<7(+l, —1, —1). The equations of motion corresponding to (2.1) are, 

D^rrza + {D,z b y{D»z b )z a = o, (2.2) 

together with Z* a Z a = 1. In this paper, we will concentrate on solutions with 
winding number k, defined by 

k = ir 1 2 e jk (DiZ a )*(D k Z a ) d 2 x, (2.3) 

equal to 2. 

The initial data are calculated from the complex functions, 
Wi(0, z, z*) = \z 2 , w 2 (0, z, z*) = /3z, 

d oWl (0, z, z*) = Xv, d Q w 2 {0, z, z*) = 0, (2.4) 

where A, (3 and v are real positive constants, and where Wi(t,z,z*) and Z a 
are related through 

i ( M 

==_ Wl . (2.5) 

+ kii 2 + |w 2 | 2 \ w 2 y 

For these initial data the energy is finite and the winding number is 2. The 
initial data describe a head-on collision of solitons which 'coincide' at t — 0. 
Negative time is the time before, and positive time is the time after the 
collision. The condition, |wi| 2 + |u> 2 | 2 = at z = 0, which holds for the 
initial data, reflects the fact that the solitons merge at t — 0. 

For doWi at t — 0, we have taken zero modes which give a head-on collision. 
This means that we don't need any extra potential energy to go through the 



3 



'ring' which forms at t = 0. For small v, there is hardly any excess energy, 
and the slow motion approximation, or collective coordinate approximation, 
can be used. The idea of this approximation [DJ is to describe the solution 
at each fixed time in terms of a configuration with minimal potential energy. 
Then the action is minimized to obtain the collective coordinates, which 
parameterize these minimal-energy configurations, as functions of time. 

For the CP 2 model the configurations with minimal potential energy, 
which for constant parameters are time independent solutions, are of the form 
(2.5) where w± and w<i are rational functions of z |14|- These configurations 
can also be given in the form, 

Z a = ^ with \f\ = 4W a , (2-6) 

where the / a 's are polynomials in z. To relate the two descriptions use is 
made of the U(l) symmetry of the theory. In terms of f a , the equations of 
motion (2.2) read 

i/i 2 " i/i 2 m 4 ■ 1 j 

The winding number can be expressed as 

1 r 2n f*zd f a 
£;=--/ lim JaZ °;J d9, (2.8) 

27TJ0 r-+oo |/|2 ' v > 

where z = re 10 . This formula shows that for k = 2, we can restrict our 
attention to quadratic polynomials /„. 

In sections 4 and 5, we will also study the effect the generalized Hopf 
term, 

C Hop f = e>»'P(D lt Z a )*(D 1 ,Z a )(Z l "d p Z b ), (2.9) 

has on the scattering process. For the Lagrangian £ + nCnopf, the equations 
of motion read, 

(5 a b - Z a Z b *)D^Z h - 2Ke^ p D p Z a ((D^Z b )*D u Z b ) = 0, (2.10) 
where, as before, Z*Z a = 1. 
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3 The collective coordinate ansatz 



For the collective coordinate approximation, we take the static solutions, 
Z a (z, c r ), parameterized by parameters c r , make the parameters time depen- 
dent and minimize within this ansatz. If we work with configurations of the 
form (2.6), this means that 

n 

fa(t,z)=Y.C n a (t)z n i (3-1) 
1=0 

where n = 2, for a 2-soliton configuration. 

Within the ansatz, the potential energy is a constant and the CP 2 kinetic 
energy is 

T = (D°Z a )*(D Z a ) 

1 ((d f:)(d o n--^f:(d nf b (d°r b )). (3.2) 



This yields, 
with, 



l/l 

| r2 T d 2 x = A^cTc n b = L (3.3) 



AT = / r2 j^(5 ab z* m z n - ±- 2 c:c s b *z n+r z* m+s )d 2 x. (3.4) 

Before we minimize the action, Pfj^ Ldt, in the next two sections, we will 
simplify the ansatz (3.1). Our next arguments work for all initial data, in- 
cluding ours in (2.4), which satisfy the following conditions: f a = f a (z),a = 
1,2,3, with d z f a + for a = 2,3; d a f k = 0,k = 1,3; d f 2 = d f 2 (z) ^ 0. 
For k = 0, we find f 1 d 2 f 3 — pdlf 1 — at t — 0. To derive this result, eq. 
(2.7), the conditions on the initial data and the formulas — did 1 = d z d z * and 
-2(9i/ a )97 6 = (d z f a )d z *f b + (d z *f a )d z f b for a, b = 1,2,3 are used. Now 
f 1 can be made real by a gauge transformation, and multiplying all / a 's by 
I// 1 we can achieve f 1 = 1 for all t. So at t = we have 9q/ 3 = 0. We show 
next that d^f 3 = for t = 0. For this we need again that the initial data 
depend on z only. 

To go further we need next that d 2 f 2 at t — depends on z only. This is 
not the case for the exact solution. In fact we obtain 

d 2 of = ^rj;(d f 2 ) 2 (3.5) 
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at t — 0. However, within the confines of the collective coordinate approxi- 
mation, we are neglecting contributions due to a z* dependence of f a , and in 
that case we can actually show that all time derivatives of f 3 vanish at t — 0. 
The reasonable assumption of analyticity now leads to the time independence 
of f 3 for t > 0. 

So far we have f 1 = 1, / 3 = (3z and 

, 2 = ai(t)z 2 + a 2 (t)z + a 3 (t) 

a 4 (t)z + a 5 (t) 1 ' ' 

If 04 depends on time the kinetic energy, / R2 Tod 2 x, diverges. (The inter- 
pretation of this divergence is that a change in 04 costs infinite energy.) So 
a 4 must be constant, and therefore zero because of the initial data. The 
kinetic energy also diverges for time dependent a 1 . Hence we have a\ — A 
and f 2 = \z 2 + 5(t)z + ^(t). Next we want to show that 5 can be taken to 
be zero. 

The Euler-Lagrange equations for L defined in eq. (3.3), the geodesic 
equations, are, 

a A mn a A mn a A rn 

Amn--n , Ujtl ab -r-n , Ujtl ab -r*-n _ u - rL cb -r*-n _ n C\7\ 

ab b dc r c dc c* dc T ~ 

with c™(t) defined in (3.1). If c°(t) = 7(t) and the other c™ are constant, then 
the equations reduce to, 

fi 4 m0 

A^l + = 0. (3.8) 



Here we have used that, 

dc 0* Q c m* ■ 

That (3.9) holds can be seen by taking the corresponding derivatives of the 
integrands which define the A's in eq. (3.4). 

If on the other hand we first set c% = ^(t) and the other c™ constant in 
eq. (3.3), and then minimize, we obtain, 

a /too 

< 7 + ^7 2 = 0. (3.10) 
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The equations (8) are now satisfied if, 

r)A 00 r)A m0 

In our case, so far the ansatz has been reduced to f 1 — 1, f 2 — \z 2 + 5(t)z + 
j(t), f 3 = f3z. For this ansatz, we have to solve the geodesic equations (3.7). 
A special solution is provided by 5 = 0, which satisfies the initial condition, 
and a function 7 which satisfies (3.10). Now (3.11) holds trivially, since 
both the integrand of and the integrand of dA^/d^f have the symmetry 
1(9 + it) = —1(9), as functions of the angle 9. Hence the corresponding 
integrals vanish. Therefore, finally our collective coordinate ansatz is 

f = l, f = \z 2 + 1 (t), f = (3z. (3.12) 

If the Hopf term is included, the Lagrange function, defined in (3.3), will 
aquire the extra term k J r2 £Ho P fd 2 x, where 

C Hopf = e^(d,Z* a )(d u Z a )Z* b d p Z b 

= ^mraMrMdof - f%n) (3.13) 

We will study the effect of the Hopf term for k small enough so that the 
ansatz (3.12) still gives a reasonable approximation. If k is not small, we do 
not expect a time independent fs to be a good approximation because 

2 ^(3v\X\ 2 z*(2+\(3\ 2 \z\ 2 ) 

doh ~ (TTWW+WWT 2 { ] 

holds for our initial data at t — 0. 



4 Approximate solution for solitons close to- 
gether 

The ansatz (3.12) leads to 

1 \\z 2 + 'y\ 2 
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Note that T is a rational function of r, where z = re id , and that the denom- 
inator is a product of quadratic functions in r 2 . So the r integration can be 
done, leading to functions of 6 which have to be integrated numerically. We 
also obtain, 

C Hop f = ~ A*^ 2 )7 - (7 - (4.2) 

Again this is a rational function of r which can be integrated over r. 

In this section, we concentrate on solitons close together, i.e., we choose 
|7| to be small. For small I7I, we expand, 

\yz 2 + a*72* 2 +0(H ^ 



|/| 2 1 + |/3| 2 r 2 + |A| 2 r 4 (1 + |/3| 2 r 2 + |A| 2 r 4 ) 2 
Now we obtain, 

/ T d 2 x = A^* + O(\ 1 \ 3 ), (4.4) 

where 

r(l + |/3| 2 r 2 ^ 
/o (1 + |/3| 2 r 2 + |A| 2 r 4 ) 
The integral can be evaluated and given in terms of 



9(A,/J) = 



This yields, 



^M^ + «) /or 2|A|< |/3 | 2 

T*f for 2\\\ = \P\ 2 (4-6) 
^4 -(f-arctan , % - ) /or 2|A| > |/3| 2 



^ = 4|A|2 ^ |/3|4 [|/^| 2 - (|/3| 4 - 2|A| 2 )^(A,^)] (4.7) 



for |/5| 2 + 2|A|. For |/3| 2 = 2|A|, the formula reduces to A = 2tt/(3|A|). 
We also have, 



/ C Hopf d 2 x = .5(7*7 - 77*) + 0(| 7 | 3 ), (4.8) 



where 



B = *f Jr ( : + ffif 1 (49) 



Note that (3 = (CP 1 embedding) implies 5 = 0. Again the integral can be 
evaluated, and we obtain 

B = m ^l w? M\W - 7|A| 2 ) + |A| 2 (10|A| 2 - W)g{\P)\ (4.10) 

for \/3\ 2 ^ 2|A|. If \(3\ 2 = 2 j Aj , the result reduces to B = llvr/60. 
For small |7|, we end up with the Lagrange function, 

L = A77* +«5(f7 -77*). (4.11) 

Expressing 7(t) = R(t)e la ^ in terms of real functions, we obtain, 

L = A(R 2 + R 2 & 2 ) — 2kBR 2 &, (4.12) 

and the equations of motion, 

AR — ARb 2 + IkBRcf = 0, (4.13) 

AR 2 '& + 2ARR& - 2kBRR = 0. (4.14) 
The second equation implies, 

ci + kBR 2 

If ci 7^ 0, 6" — > 00 as i? — > 0. So we set c\ = 0, and obtain & = kB/A and 
w.l.g. a = nBt/A. Now the equation for R reads, 

A 2 R + k 2 B 2 R = 0, (4.16) 

and, using R(0) = 0, we obtain, 

„ KBt kB 

R = c 2 sm— — = c 2 —^-t, (4.17) 

for small R (and t). 

We now consider |/|. At t — 0, |/| is radially symmetric and describes 
a 'crater-like' surface. When t becomes nonzero, the surface 'buckles'. Con- 
sider, with A real for simplicity, 

|/| 2 = 1 + |/5| 2 r 2 + |A| V + 2Xr 2 R(t) cos(2fl - a(t)). (4.18) 
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The 9 derivatives gives us the 'ridges' and 'valleys', which are at 6 = nn/2 + 
a/2. If for definiteness A > and c 2 > 0, the 'valleys' are at ^ and 7r+ ^ 
for t < 0, and at f + ^ and ^ + ^ for t > shortly before and shortly 
after the collision. For k — 0, we again find 90° scattering, whereas for k^O, 
we find a deviation from this type of scattering. 



5 Solitons at a distance 



To study solitons which are not close together, T , from Eq.(4.1), and Cnopf, 
from Eq.(4.2), must be integrated over R 2 . For T , the result of the r inte- 
gration can be written as 

r /*oo 

/ T d 2 x= / d0a(A,/3,7,0)77*. (5.1) 
Jr. 2 Jo 

a(A, /?, 7, 0) can be expressed in terms of 

6(A, P, 7 , 0) = |/3| 2 + \l*e™ + A* 7 e- 2 * e , (5.2) 
A(A,/3, 7 ,^ = 4|A| 2 (l + | 7 | 2 )-6 2 , (5.3) 

and 

7=s l0 §£v9 f° r A<0 
h(\,[3, 1 ,9) = { f /or A = (5.4) 

-^(f-arctan^) for A>0 



as 

m „ m l/3| 2 2|A| 2 -|/5| 2 6 

a(A '^'^ g ^ 4|A| 2 (l + |7| 2 ) + ^AFA^ 



2|A| 2 MA,/5,7,0)- 



1 + 



(5-5) 

The integration and the integration of the Euler-Lagrange equations have 
to be done numerically. 

For large |7| though, we can obtain analytical results. We start by inte- 
grating the factor of 77* on the right-hand side of (4.1) over regions close to 
±a. Set 7 = — Aa 2 and z = ±a + b with < |6| < |a| 3 / 4 . (Here and in our 
arguments below, there is a wide choice of limits of integration. Here for 
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any power of \a\ between | and 1 would do). Neglecting b terms relative to 
a terms, we get the integrals, 

■ r 2 i+ i^i 2 h 2 m * (5 6 ) 

Jo (1 + \3\ 2 \a\ 2 + MW 2 \a\ 2 \b\ 2 ) 2 ' ' ' 1 ^ 



47T 



Hence, 

I± = 

and for large \a\, 



7T 



tt(1+ |/3| 2 |a| 



|A| 2 |a| 2 |A| 2 |a| 2 (l + |/5| 2 |a| 2 + 4|A| 2 |a| 7 / 2 )' 



(5.7) 



^^ +I -ww- (5 - 8) 

We now show that the rest of the integral is of lower order than (5.8) and 
can be neglected. This is an integral over R 2 with two circles centered at ±a 
taken out. We split the integral into three integrals, I±, I2 and 1%. I\ is the 
integral for \z\ between and |a| 3//4 . I2 results from the integration for \z\ from 
|a| 3 / 4 to |a| 2 , and I3 is the integral over the rest of R 2 outside the big circle 
of radius \a\ 2 . Starting with the integral I±, we use the inequalities \z ± a\ > 
\a\/2, and take the lower bounds of these inequalities in the integrand. This 
yields, 

T < 4 f H3/2 1 + IH 2 M 2 ,2 

1 " 71 Jo (l + |/5| 2 |z| 2 + |A| 2 |a| 4 /16) 2 lZl 

So the bound on I x is of order l/|a| 5 , and I\ can be neglected. 

The second integration area, which is an annulus with two small disks 
deleted, is split into two symmetric halves with a small disk deleted in each. 
In one half, \z — a\ > |a| 3 / 4 and \z + a\ > \z\; in the other, \z + a\ > |a| 3//4 
and \z — a\ > \z\. By taking the lower bounds of these inequalities in the 
integrand, we get an upper bound for I 2 . We now increase the value of this 
upper bound by including the two small disks. So we have 

H 4 1 + |/3| 2 M 2 „ l2 4tt 



x 



/2 - 47r i,a/ 2 (1 + |^| 2 N 2 + |A| 2 - — dlzl = 
|/5| 2 log(l + |/5| 2 |z| 2 + |A| 2 |a| 3 / 2 



|a|3/2 (1 + \[3\ 2 \z\ 2 + |A| 2 |a| 3 / 2 |^| 2 ) 2 1 1 (\(3\ 2 + |A| 2 |a| 3 / 2 ) 2 



\M 2 


a| 3 / 2 


i + \P\ 2 


z 


2 +|A| 2 


a\ 3 / 2 \z\ 2 



J |a| 3 /2 

(5.10) 
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For large \a\, this goes like , and J 2 can be neglected. Finally, for the 
integral J 3 , we use the inequalities |z±a| > \z\/2. This leads to the following 



upper bound, 
h < 4vr 



l-l 4 (1 + 



+ |A| 2 |^| 4 /16) 



-Az\ 2 < 



2 w n\P\ : 
\X\ 4 \a\ 8 



(5.11) 



To include the generalized Hopf term, we again perform the r integration 
and obtain 



Jn 2 



d 2 xC 



Hopf - 



I 2 * d9c(X, 0,7,7,0) , (5-12) 
Jo 



where 



c(A,/3, 7 ,7,0) 



7*7 - 77* + 



21 A| 



(A* 7 e- 2ie - A 7 *e' 



X 



2IAI 



■h(\,P,>y,0)- 



(5.13) 



A v " ' A(l + | 7 | 2 )_ 

Again, the 9 integration and the integration of the Euler-Lagrange equations 
have to be done numerically. 

For large \a\, however, we obtain some analytical results. We must in- 
tegrate Cnopf, given in (4.2), over R 2 . For large \a\, the leading term will 
again come from the integration over regions close to ±a. As above we set 
z = ±a + b and integrate for < \b\ 2 < |a| 3 / 2 . Neglecting b terms relative to 
a terms, these integrals can be written as, 



J. 



± 



Jo 



4m\P\ 2 \\\ 


2 \a 


2 ( 


a* 


a 




ad 




(1 + I0I 2 


a\ 


2 + 4|A 


2 


a 


2 


6 |2)2 



d\b[ 



To leading order in \a\, we obtain, 



J=J + + J_ = 



2m 



(a*d — ad*). 



(5.14) 



(5.15) 



We now show that the integral over R 2 , with the two disks centered at ±a 
taken out, can be neglected compared to J. Again using the circles of radius 
\z\ = |a| 3 / 4 and \z\ = \a\ 2 , respectively, the integral is divided into three 
integrals, Ji, J 2 and J 3 . With the help of the inequalities \z ± a\ > \a\/2 
which we used above, we now get the bound, 



l^il < 



2 




d 




A 


2 


a 


7/2 



(5.16) 



12 



on the integral over the disk at the origin. J\ is clearly of lower order than 
J. Integration over the annulus for | cz | 3 / 2 < \z\ 2 < |a| 4 , with the two disks 
deleted, leads to the bound, 

\U < 



4 7 r| / 3| 2 |A| 2 |a| 


a 




m 2 +\M 2 


|a| 3 / 2 ) 2 



M 4 

log(l + |/?| 2 |d 2 + |A| 2 |a| 3 / 2 |d 2 ) + 



(1- 


a 


w- 




A 


2 


a 


7/2) 


1 + 1/3 


I 2 


z 


2 + 


A 


2 


a 


3/2 |;j 


2 



X 

! ' ' I 1 2 I - 1 2 i I \ 1 2 I , , I : i / 2 1 - 1 2 

|a| 3 / 2 

(5.17) 

This bound is of order |a|/|a| 3//2 , and J 2 can be neglected. J 3 can be neglected 
as well, since 

, , 2 10 7r|/?| 2 |a| , s 

If we use I7I 2 = 4|A| 2 |a| 2 |d| 2 and eqs. (5.8) and (5.15), we find the follow- 
ing Lagrange function for large \a\, 

L = 8naa* + 2mn (- - — V (5.19) 

\a a* J 

The Hopf term is the time derivative of 2itik In ^- and does not contribute 
to the equations of motion. To leading order in the distance between the 
solitons, we have a(t) = a Q + vt. This is the free motion seen in previous 
numerical studies. 



6 Conclusions 



In this paper, which a sequel to |TT|, we continued our investigation of the 
scattering of solitons in the CP 2 model in (2+1) dimensions. In [ |TT[ 



we 



dimensions. In 

showed that the addition of an extra term, the "generalised Hopf term" 
changes dynamics of the solitons; the 90° scattering seen in head-on colli- 
sions of CP 2 solitons in the absence of this term is replaced by a difnection, 
which for small values of the coefficient of this term is proportional to this 
coefficient. These expectations, based on the expansion of the field at the 
time of overlap of the two solitons, were shown to be be born out by the 
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results of numerical simulations. The simulations have also shown that at 
larger distances the generalised Hopf term seems to play little role; all de- 
flection of the solitons seems to be confined to the region when they are very 
close together. 

In this paper we have looked at this problem from the point of view 
of collective coordinates. We have argued that in this problem we can ap- 
proximate the field configuration by a simple polynomial expressions (with 
many terms vanishing) and then looked at the effects of the nonvanishing 
coefficient of the generalised Hopf term on the dynamics of these collective 
coordinates. We have shown that this approximation is self consistent and 
that it captures the main features of the data (ie the deflection from the 90° 
scattering beeing proportional to the coefficient of the Hopf term). Supris- 
ingly, the approximation also captures the effect of the Hopf term becoming 
irrelevant (for classical dynamics) when the solitons are well separated; in 
this case its contribution to the effective lagrangian becomes a total deriva- 
tive. Moreover, the form of this total derivative is very similar to what is 
seen in the CP 1 model (though there this form is valid for all distances be- 
tween solitons, including the cases when they overlap). Thus we see that well 
separated solitons in CP 2 model are not that different from solitons of the 
CP 1 model; only when they are close together the nature of target manifold 
begins to play an important role, and the generalised Hopf term brings out 
this difference. 

In addition our work has given an extra support for the use of collective 
coordinates; they provide a good approximation to the full dynamics also 
when the generalised Hopf term is present. 
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